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Abstract: We consider two non-standard cosmological scenaria according to which the 
universe is reheated to a low reheating temperature after the late decay of a scalar 
field or is dominated by the kinetic energy of a quintessence field in the context of 
a tracking quintessential model. In both cases, we calculate the relic density of the 
Weakly Interacting Massive Particles (WIMPs) and show that it can be enhanced with 
respect to its value in the standard cosmology. By adjusting the parameters of the low 
reheating or the quintessential scenario, the cold dark matter abundance in the universe 
can become compatible with large values for the annihilation cross section times the 
velocity of the WIMPs. Using these values and assuming that the WIMPs annihilate 
predominantly to yU + //~ or t + t~, we calculate the induced fluxes of e ± cosmic rays 

and fit the current data of PAMELA and ATIC or Fermi LAT. We achieve rather good fits 
especially to PAMELA and Fermi-LAT data in conjunction with a marginal fulfillment 
of the restriction arising from cosmic microwave background, provided that the WIMPs 
annihilate predominantly to ju + yu~ '. In both non-standard scenaria the required transition 
temperature to the conventional radiation dominated era turns out to be lower than about 
0.7 GeV. In the case of the low reheating, an appreciable non-thermal contribution to the 
WIMP relic density is also necessitated. 
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1. Introduction 

The accurate determination of cosmological parameters by up-to-date observations, most notably 
by the Wilkinson Microwave Anisotropy Probe (WMAP) [1,2], establishes a quite extensive and 
convincing evidence for the constitution of the present universe by an enigmatic component called 
Cold Dark Matter (CDM) with abundance, Qcdm^ 2 , in the following range: 

^cdm/i 2 = 0.1099 ± 0.0124 (1.1) 

at 95% confidence level (c.L). Natural candidates [3] to account for the CDM are [4] the weakly 
interacting massive particles (WIMPs, hereafter denoted as x) with prominent representative (for 
other WIMPs, see Ref. [5, 6]) the lightest neutralino [7] which turns out to be the lightest super- 
symmetric (SUSY) particle (LSP) in a sizeable fraction of the parameter space of the SUSY models 
and therefore, stable under the assumption of the /^-parity conservation. In view of Eq. (1.1), the 
relic density of ^'s, Q. x h 2 , is to satisfy a very narrow range of values: 

(a) 0.097 < Q. x h 2 and (b) Q x h 2 < 0.12, (1.2) 

with the lower bound being valid under the assumption that CDM is entirely composed by ^f's. 

The calculation of Q. x h 2 crucially depends [8] on the adopted assumption about the dominant 
component of the universe during the decoupling of WIMPs. The usual assumption is that this 
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occurs during the radiation dominated (RD) epoch which commences after the primordial inflation. 
However, our ignorance about the universal history before Big Bang nucleosynthesis (BBN) allows 
for other possibilities. E.g., the presence of a scalar field, which dominates the budget of the 
universal energy density through its potential [8, 9, 10, 11, 12, 13] or kinetic [8, 14, 15, 16] energy 
density, can enhance significantly Q. x h 2 with respect to (w.r.t) its value in the standard cosmology 
(SC). In the first case, the scalar field can generate an episode of low reheating which can be 
accompanied by thermal and/or non-thermal production of %'s. In the second case, a kination 
dominated (KD) epoch [17], which may be embedded [16, 18, 19, 20, 21] in a quintessential 
framework, can arise. As a bonus, in the latter case, the problem of the second major component of 
the present universe, called Dark Energy (DE) can be addressed - for reviews see, e.g., Ref. [22]. 

The aforementioned enhancements of Q. x h 2 have attracted much attention recently [23] since 
they assist us to interpret, through WIMP annihilation in the galaxy and consistently with Eq. (1.1), 
the reported [24, 25] excess on the positron (e + ) and/or electron (e~) cosmic-ray (CR) flux, with- 
out invoking any pole effect [26], ad-hoc boost factor [27] or other astrophysical sources [28]. In 
particular, PAMELA experiment has reported [24] (confirming previous experiments [29]) an un- 
expected rise of e + flux fraction for values of the e + energy, E e +, in the range (10 - 100) GeV, in 
contrast to the power-law falling background. Moreover, data by the ATIC experiment [25] shows 
an excess in the total e + and e~ flux for 300 < E e + /GeV < 800. On the other hand, the very recently 
released data from Fermi LAT indicates [30] smaller fluxes than the ATIC data in the same range 
of energies. Nevertheless, we consider (separately) both latter data in our study. 

In this paper we reconsider the increase of Q. x h 2 within a low reheating scenario (LRS) or 
a quintessential kination scenario (QKS) in light of the experimental results above. Namely, we 
recall (Sec. 2) comparatively the salient features of the two non-standard scenaria, solving nu- 
merically the relevant system of equations, reviewing the cosmological dynamics and imposing 
a number of observational constraints. Particularly, in the LRS we consider the late decay of a 
massive field which reheats the universe to a low reheating temperature. In the QKS, we consider 
the recently implemented [21] generation of a KD era (associated with an oscillatory evolution 
of the quintessence field) in the context of tracking quintessential model with a Hubble-induced 
mass term for the quintessence field. We then (Sec. 3) investigate the enhancement of Q. x h 2 w.r.t 
its value in SC within these non-standard scenaria. We show that the increase of Q. x h 2 depends 
on (i) the reheat temperature and the number of ^'s produced per decay and unit mass of the de- 
caying field, in the case of LRS, and (ii) the proximity between the freeze-out temperature and the 
temperature where the evolution of the quintessence develops extrema, in the case of QKS. We 
also present (Sec. 4) the energy spectra of the e ± -CR, assuming mate's annihilate into e + e~, 
or t + t~ and adopting an isothermal halo profile [31]. Although Cosmic Microwave Background 
(CMB) [32, 33, 34] tightly constrain the relevant parameters, we achieve rather satisfactory fittings 
especially to the combination of PAMELA and Fermi-LAT e ± -CR data and for the case where ^'s 
annihilate into jU + //~ (Sec. 5). Fulfilment of Eq. (1.2) is also possible by appropriately adjusting the 
parameters of the LRS or QKS. We end up with our conclusions in Sec. 6. 

Throughout the text, brackets are used by applying disjunctive correspondence, the subscript 
or superscript is referred to present-day values (except for the coefficient Vo) and log [In] stands 
for logarithm with basis 10 [e\. Besides Sec. 4, natural units for the Planck's constant, Boltzmann's 
constant and the velocity of light (h = c = &b = 1) are assumed. 
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2. Non-Standard Cosmological Scenaria 

In this section we present comparatively the main features of the two non-standard cosmological 
scenaria under consideration. Namely, in Sec. 2.1 we expose the basic assumptions of each scenario 
with reference to the SC and introduce notation. Sec. 2.2 and Sec. 2.3 are devoted to a review of the 
LRS and QKS respectively. Despite the fact that the displayed scenaria have been already analyzed 
in Ref. [11, 21] we prefer to briefly recall and update our results for completeness and clarity. 

2.1 The General Set-up 

According to SC, primordial inflation is followed by a RD era. The x species (i) are produced 
through thermal scatterings in the plasma, (ii) reach chemical equilibrium with plasma and (iii) 
decouple from the cosmic fluid at a temperature Tp ~ (10 - 20) GeV during the RD era. The 
assumptions above fix the form of the relevant Boltzmann equation, the required strength of the^ 
interactions for thermal production (TP) and lead to an isentropic cosmological evolution during 
the^f decoupling: The Hubble parameter is H oc T 2 with temperature T oc R 1 where R is the 
scale factor of the universe. In this context, the Q. x h 2 calculation depends only on two parameters: 
The,Y mass, m x and the thermal-averaged cross section of x times velocity, (cry) - see Table 1. 
Although (crv) can be derived from m x and the residual (s)particle spectrum once a low energy 
theory is adopted (see, e.g., Ref. [15]), we treat m x and (crv) as unrelated input parameters in order 
to keep our approach as general as possible (see, e.g., Ref. [16, 31]). Also, to be in harmony with 
the assumptions employed in the derivation of the restrictions mentioned in Sec. 5.1, we consider 
throughout constant (or>)'s, i.e., independent of T. 

The modern cosmo-particles theories, however, are abundant in scalar massive particles (mod- 
uli) which can decay out of equilibrium when H becomes equal to their mass creating episodes of 
reheating. In the LRS, we assume that such a scalar particle <f>, with mass m^, decays with a rate 
into radiation, producing an average number N x of ^'s. The key point in this case is that the re- 
heating process is not instantaneous [9, 11]. During its realization, the maximal temperature, T max , 
is much larger than the so-called reheat temperature, Trh, which can be taken to be lower than Tp. 
Also, for T > Trh, H x T 4 with T x 7?~ 3 ^ 8 and an entropy production occurs (in contrast with 
the SC). The x species (i') decouple during the decaying-^ dominated era (ii') do or do not reach 
chemical equilibrium with the thermal bath (iii') are produced by thermal scatterings and directly 
by the <f> decay (which naturally arises even without direct coupling [9]). As a consequence, the 
Q. x h 2 calculation depends also on Trh, and N x - see Table 1. 

Another role that a scalar field could play when it does not couple to matter (contrary to <f>) is 
this of quintessence. In our QKS, such a scalar field q (not to be confused with the deceleration 
parameter [2]) is supposed to roll down its inverse power-law potential (with exponent a and a mass 
scale M) supplemented with a Hubble-induced mass term (with coefficient b) motivated mainly by 
non-canonical Kahler potential [20, 35] - c.f. Ref. [36]. A mild tuning of b and of the initial 
conditions at an initial temperature T\ - Hi = H{T\) and qi = q{T{) - ensures the coexistence of 
an early modified KD phase with the tracking [37, 38] solutions and offers the desirable property 
of the insensitivity to the initial conditions [20, 21]. Since the q kinetic energy, which decreases as 
T 6 (except for isolated points) dominates we get H x T 3 with T oc R l . If the ^-decoupling occurs 
during this KD phase - the assumptions (i) and (ii) are maintained - the Q. x h 2 calculation depends 
also on a, b, M, Hi, Ti and qi in this scenario - see Table 1. 
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Table 1 : Comparing the SC with the LRS and the QKS (the various symbols are explained in Sec. 2.1, the subscript I is 
referred to the onset of each scenario and "est" stands for "constant"). 



In the two non-standard scenaria under consideration, H is given by 

H 2 = J_((p<p+P x + Pr) for the LRS, ^ ^ 



3m 2 \ (p q + p R + p M ) for the QKS, 

where p, with i - <p,q andx is the energy density of <p, q and^ respectively and mp = Mp/ V8tt 
where M P = 1.22 • 10 19 GeV is the Planck mass. The energy density of radiation, pr, and the 
entropy density, s, can be evaluated as a function of T, whilst the energy density of matter, pm, 
with reference to its present-day value: 
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PR = ^Q#p* r4 ' s = -^-9s*T 3 and PmR 3 =Pmo#o (2.2) 

where q p *(T) [g s *(T)] is the energy [entropy] effective number of degrees of freedom at temperature 
T. Their precise numerical values are evaluated by using the tables included in public packages 
[39] and assuming the particle spectrum of the Minimal Supersymmetric Standard Model. 

The initial value of H, Hi, in both non-standard scenaria can be restricted, assuming that a 
primordial phase of inflation (driven by a scalar field different from (f> or q, in general) is responsi- 
ble for the generation of the power spectrum of the curvature scalar P s and tensor P t perturbations. 
Indeed, imposing the conservative restriction r = P t /P s < 1 and using the observational [1] nor- 
malization of P s , an upper bound on H\ can be found as follows: 



n 



Hi < -—m P P l J 2 =* Hi% 2.65 • 10 14 GeV, (2.3) 
V2 

where * means that P s * is measured at the pivot scale k* - 0.002/Mpc. 

Let us, finally, introduce a set of normalized quantities which simplify significantly the relevant 
formulas. In particular we define 

Pi - Pi/PcQ, with i = R, M, cf> and q, J - J/p 3 C Q with J = n x , n x q and n$, (2.4a) 
rhj - mil p l J£ with i = x and <p,J = J/Hq with J - H and and (cry) - V3mpp^ 4 (crv) (2.4b) 

where itj with i = <f> and x is the number density of x and <f> respectively. Note that p x - m x n x and 
Ptp = A^«0 where A^ - (m^ - N x m x )lm^,. In our numerical calculation, we use the values: 

p c0 « 8.1 • l(T 47 /z 2 GeV 4 , with h = 0.72, p M o - 0.26 and T = 2.35 • 10~ 13 GeV. (2.5) 

We have also H = 2.13 • 10 _42 /i GeV and from Eq. (2.2) we get p R0 - 8.04 • 10" 5 . 
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2.2 The Low Reheating Scenario 

We below (Sec. 2.2.1) present the system of equations which governs the cosmological evolution 
in the LRS, summarize (Sec. 2.2.2) the various observational restrictions that have to be imposed 
and sketch (Sec. 2.2.3) the basics of the relevant dynamics. 

2.2.1 Relevant Equations. Under the assumption that the decay products of <j> are rapidly ther- 
malized (see below) the energy densities p^ and pr obey the following Boltzmann equations 

p^ + 3flp0 + r p = 0, (2.6a) 
p R + 4Hp R - I> - 2m x {o-v) {n\ - nf) = 0, (2.6b) 

where dot stands for derivative w.r.t the cosmic time, t and n x q is the equilibrium number density of 
X, which obeys the Maxwell-Boltzmann statistics: 

is obtained by expanding the modified Bessel function of the 2nd kind of order n for x «; 1. 
Assuming that x's are Majorana fermions, we set g - 2 for their number of degrees of freedom. 
Note that although in our numerical program Eqs. (2.6a) and (2.6b) are resolved together with 
Eq. (3.2) - see Sec. 3.1-, we here prefer to present just the two first equations since the influence 
of n x to the dynamics of reheating via the last term of the left hand side in Eq. (2.6b) is in general 
negligible. Moreover, can be replaced by Trh through the relation [11]: 



r ^ 5 V^5 = V 72 W (2 " 8) 

Note that the adopted prefactor, which is slightly different than those used in the bibliography 
[10, 12], assists us to approach accurately the numerical solution of p^(rRn) = Pr(7rh). 

The numerical integration of Eqs. (2.6a) and (2.6b) is facilitated by absorbing the dilution 
terms. To this end, we find it convenient to define the following variables [10, 11]: 

U = PV? 3 , A = PrR\ fx = "x r3 and fx q = n ? R3 ( 2 - 9 ) 
and convert the time derivatives to derivatives w.r.t the logarithmic time [11]: 

f = ln(R/R{) => R' =R and R = R^ (2.10) 

where prime in this section denotes derivation w.r.t f and the value of R\ in this definition can be 
conveniently selected so as the resolution of the system is numerically stable. After realize the 
modifications above, Eqs. (2.6a) and (2.6b) become: 

(a) Hf; = -t^U and (b) HR 2 f R = T^R 3 + 2m x W)(f x - ff) , (2.1 1) 

where H and T can be expressed correspondingly, in terms of the variables in Eq. (2.9), as follows: 

/ / 30 f 

(a) H = R- 3 ' 2 J ft + f R /R and (b) T - \-^-^ Pc0 j . (2.12) 
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The system of Eq. (2.1 1) can be solved from to f f ~ 50, imposing the following initial conditions 
(recall that the subscript I is referred to quantities defined at x = 0): 

H\ = m.$ => pfa = m^/Hl and p RI = p Xl - 0. (2.13) 

However, the results on Q. x h 2 do not depend on the explicit value of p^, as long as Trh < Tp < T max , 
and are invariant [11, 12] for fixed N x m~ l (and T$, m x , (crv)). Therefore, for presentation purposes, 
it is convenient to define the following quantity [12]: 

100 TeV 

( 2 - 14 ) 

2.2.2 Imposed Requirements. We impose on the LRS the following requirements: 

• The BBN Constraint. The presence of p^ should not jeopardize the successful predictions of 
BBN which commences at about r BBN = 1 MeV [40]. Namely, we require: 

r RH > 1 MeV (95% c.L). (2.15) 

Given that tp decays mostly through gravitational interactions, T^ and consequently Trh - see 
Eq. (2.8) - are highly suppressed. Therefore [12] fulfilment of BBN constraint with more or 
less natural coupling constants requires > 100 TeV. 

• Constraints on the range of m^. Eq. (2.3) assists us to impose an upper bound on due to 
our initial condition in Eq. (2.13). On the other hand, can be bounded from below, too, 
demanding the decay of <p to a pair of x's (with mass m x ) to be kinematically allowed. All in 
all we require: 

2m x < m(P < 2.65 • 10 14 GeV. (2.16) 

Note that the upper bound of Eq. (2.16) assures also the rapid thermalization of the 0-decay 
products. Indeed, the latter condition, which is crucial for Eqs. (2.6a) and (2.6b) to be appli- 
cable, is satisfied [41] for m < 8 • 10 14 GeV. 

• Constraint on the range of N x . Depending on the coupling between <f> andx in a specific theory, 
a variety of N x s is possible [9, 12, 42]. In our approach we conservatively set the upper bound 
N x < 1. 

Let us finally mention that, on quite general ground, any modulus (f> has an unsuppressed 
coupling to gravitino, G. Possible decay of cf> to G creates the co-called moduli-induced G problem 
[42]. To avoid these complications, we are obliged to assume that the masses of G and (f> are of the 
same order of magnitude. 

2.2.3 Dynamics of Reheating. The cosmological evolution of the various quantities involved in the 
LRS as a function of x is illustrated in Fig. 1 for - 100 TeV and (A^, Trh) = (1, 0.5 GeV) (gray 
lines) or (N x , Trh) = (7.5 • 10~ 5 , 1 MeV) (light gray lines). In particular, we design log p,- with i - <p 
(solid lines) and i =R (dashed lines) [T] versus x in Fig. 2-(a) [Fig. 2-(b)]. The quantities p, with 
i = <f> and R [T] are computed by substituting the numerical solution of Eq. (2.11) into Eq. (2.9) 
[Eq. (2.12b)]. From Fig. 1 we can understand the dynamics of the universe during the two distinct 
phases [10, 11]: 



7 



CDM in non-Standard Cosmologies, PAMELA, ATIC and Fermi LAT 



.- 60 

idr 50 

O) 

O 40 

30 
20 
10 



1 1 1 1 1 ■ 1 1 1 ■ 1 1 ■ 1 1 ■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



T RH = 0.5GeV 



T RH= 1 



MeV 



i = R 



m =100TeV 



10 15 20 25 30 35 

f = ln(R / R,) 

(a) 



40 45 50 



10= 
10 4 
10 3 
> 10 2 

a> 

O 10 r 

h 1 f 

10- 1 

10 2 t 



T„ 



:0.5GeV 



T RH = 1 MeV 



10 4 



m =100TeV 



10 15 20 25 30 35 

f = ln(R / R,) 

(b) 



40 45 50 



Figure 1: The evolution of (a) logp, with i = (/> (solid lines) and R (dashed lines) and (b) T as a function of x = \n(R/Ri) 
form* = 100 TeV and (N x , T m ) = (1, 0.5 GeV) (gray lines) or (N x , T m ) = (7.5 • lO" 5 , 1 MeV) (light gray lines). In both 
cases, wetake£l x h 2 = 0.11 form x = 0.5 TeV and {crv) = 3 • 10~ 7 GeV~ 2 . 



For T » Trh, we have » pr. Consequently, inserting H - ^/p^into Eq. (2.11) we extract: 



(a) p^ - and (b) p R = -T^pf (e^ 2 - e~ 4f ) . 



(2.17) 



The function pr(t) reaches at f max - ln(1.48) - 0.39 a maximum value pR max - 0.14 y/p^T^ 
corresponding to a T - T max derived through Eq. (2.2). The completion of the reheating is 
realized at x = trh, such that: 



PrCcrh) - P^(trh) 
where a corner is observed on the curves of Fig. l-(b) 



2 2- 
trh - In -r 



-1/2 



(2.18) 



• For T < Trh, we get pr » p^, and so, H - a/pr. Plugging it into Eq. (2.11) we can obtain 
approximately the following expressions: 

P - P^rh) exp |-3(x - f RH ) - ^ (e 2(f ^ RH) - l)j and p R = Pr(trh) g-^^n). (2. 19) 
2.3 The Quintessential Kination Scenario 

We present below (Sec. 2.3.1) the equations which govern the cosmological evolution in the QKS 
and then enumerate (Sec. 2.3.2) the various observational restrictions that have to be imposed. We 
also highlight the q dynamics (Sec. 2.3.3) and describe the allowed parameter space (Sec. 2.3.4). 

2.3.1 Relevant Equations. Under the assumption that q is spatially homogeneous, it obeys the 
following Klein-Gordon equation: 



q + 3Hq + V, q = 0, where V = V a + V b with V a = 



M 



4+a 



and V b - -H 2 q 2 , (2.20) 



is the adopted potential for the field q with M a mass-scale and , q stands for derivative w.r.t q. 
In our approach Vb is present throughout the cosmological evolution of q. The induced coupling 
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between q and CDM during the matter dominated era is too suppressed to have any observational 
consequence. Nonetheless, we have checked that our results remain intact even if we switch off this 
term after the onset of the matter domination. The numerical integration of Eq. (2.20) is facilitated 
by converting the time derivatives to derivatives w.r.t the logarithmic time [16] which is defined as 
a function of the redshift z: 

x = \n(R/R ) = -ln(l + z). (2.21) 
Changing the differentiation and introducing the following quantities: 

V a = VJ Pc0 , f q = qR 3 / V^o and q = q/ V3m P , (2.22) 
Eq. (2.20) turns out to be equivalent to the system of two first-order equations: 

f q = Hq'R 3 and Hf' q IR 3 + bH 2 q + bHH q q 2 + V a , q = 0, (2.23) 

where H 2 = - _ ^ {^fjR 6 + V a +p R + p M ) , (2.24) 

prime in this section denotes derivative w.r.t. x and M can be found from the dimensionless quan- 
tities as follows: 

hi- \a - \l/(4+a) 

M = (( V3m P ) V p c0 ) with V a = V /q a . (2.25) 

Eq. (2.23) can be resolved numerically if two initial conditions are specified at an initial x, ti 
corresponding to a temperature T\, which is defined as the maximal T after the end of primordial 
inflation, assuming instantaneous reheating. We take q(ti) = 10" 2 throughout our investigation, 
without any loss of generality (see below) and let as free parameter Hi - H(Ti). 

2.3.2 Imposed Requirements. Our QKS can be [21] consistent with the following restrictions: 

• Constraints on the initial conditions. We focus on the initial conditions which assure a com- 
plete initial domination of kination consistently with Eq. (2.3), i.e., 

(a) a\ = n q (T{) = 1 and (b) Bj £ 1.72 ■ 10 56 (2.26) 
with Q. q = p q l(p q + p R + p M ) where p q = q/2 + V. (2.27) 

• The BBN Constraint. The presence of p q has not to spoil the successful predictions of BBN 
which commences at about xbbn = -22.5 (7bbn = 1 MeV) [40]. Namely, we require: 

Qf N - n,(TBBN) < 0.21 (95%c.l.) (2.28) 
where 0.21 corresponds to additional effective neutrinos species 8N V < 1.6 [40]. 

• DE-Density and Coincidence Constraint. These two constraints can be addressed if (i) the 
present value of p q , p q o, is compatible with the abundance of DE in the universe [1] and (ii) p q 
has already reached the tracking behavior. In other words we have to demand [20] 

(a) Q q() - pgo - 0.74 and (b) d 2 V{x = 0)/dq 2 - H 2 , (2.29) 

where we restrict ourselves to the central experimental value of Cl q o, since, this choice does 
not affect crucially our results on the CDM abundance. 

• Acceleration Constraint. Any successful quintessential model has to account for the present- 
day acceleration of the universe, i.e., [1] (see also Ref. [43]) 

-1.12 < w q (0) < -0.86 (95% c.l.) with w q = (q 2 /2 - V)/(q 2 /2 + V). (2.30) 
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Figure 2: The evolution of (a) logp, with i = q (gray [dark gray] line), R+M (light gray line) and A (thin line); and 
(b) q (gray [dark gray] line) as a function ofx = ln(R/R ) forq Y = 0.01, a = 0.5, b = 0.2, T Y = 10 9 GeV (%i =s -51.16) 
M = 4.8 eV and Hi = 4.7 • 10 51 (Q% m = 0.0001) fffi = 2 • 10 53 (Q% m = 0.21)]. 



2.3.3 Quintessential Dynamics. The cosmological evolution of the various quantities involved in 
our QKS as a function of x is illustrated in Fig. 2 for q\ - 0.01, a = 0.5, b - 0.2, T\ - 10 9 GeV 
(ti a -51.2) M - 4.8 eV and H Y = 2 ■ 10 53 (Q™ N = 0.21 dark gray lines) or i^i = 2 • 10 53 
^bbn _ 0.0001, gray lines). Particularly, in Fig. 2-(a) [Fig. 2-(b)] we draw logp,- [q] versus t. For 
i - q (gray and dark gray lines), p q is computed by inserting in the last equation of Eq. (2.27) the 
numerical solution of Eq. (2.23). For i = R + M (light gray line), we show Pr + m - Pr + Pm given 
by Eq. (2.2). For i = A (thin black line), Pa is the dimensionless energy density of the attractor 
solution (see below). 

From Fig. 2 we can conclude that q undergoes four phases during its cosmological evolution 
[20,21,44]: 

• The kinetic-energy dominated phase during which the evolution of both the universe and q is 
dominated by q/2 » V. Therefore H - Hq' / ^2 - bq 2 and integrating it we obtain [21] 




(2.31) 



Obviously, for b > 0, q is set in harmonic oscillations during the KD era. In particular, q 
develops extrema at 

x ext - (2k + 1) + Ti with k = 0, 1, 2, ... (2.32) 

On the other hand, q - Hq' almost vanish for x = x ext . Therefore, at x - x ext , p R dominates 
instantaneously over q/2. As a consequence, the q oscillations become anharmonic. This 
phase terminates for x = xkr where p q = pr. For the inputs of Fig. 2 we get xrr = -25.6 
[x KR = -28.2] (or T K r = 0.02 GeV [T KR - 0.21 GeV]) for Hi = 2 ■ 10 53 [Hi = 4.7 • 10 51 ]- We 
observe that the lowest Trr corresponds to the largest Hi (and f2^ BN ). From Fig. 2-(b) we also 
remark that the height of the fifth peak of q decreases with Hi. In fact, for Hi < 4.7 • 10 51 we 
take qo < and so, q can not serve as quintessence (see below). 
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Figure 3: Allowed (gray shaded) region by Eqs. (2.26) - (2.30) in the (a) b-H x and (b) b - n™ N plane for a = 0.5, qi = 
0.01 and r, = 10 9 GeV. 

• The frozen-field dominated phase, where the universe becomes RD and p q is dominated ini- 
tially by q/2 and subsequently by V and q is stabilized to a constant value - see Fig. 2-(b). 

• The attractor dominated phase, where p q - V and pm dominates the evolution of the uni- 
verse. The system in Eq. (2.20) admits [44] a tracking solution since the energy density of the 



tracks pu until x - tm where the tracking regime terminates and pu - pa- For both Hi's used 
in Fig. 2, we get XAf = -0.4 whereas the onset of this phase occurs at XAi = -4.8 [tAi = -8.6] 
for Hi = 2 ■ 10 53 [Hi = 4.7 • 10 51 ]- We observe that although the used #i's differ by two 
orders of magnitude, both p q 's reach Pa highlighting thereby the insensitivity of our QKS to 
the initial conditions. 

• Vacuum Dominated Phase. For x > Tm, the evolution of the universe is dominated by V. For 
the parameters used in Fig. 2 we get w q (0) -0.88 and £l q (0) - 0.74. 

2.3.4 Allowed Parameter Space. The free parameters of our QKS listed in Table 1 can be re- 
stricted using the criteria presented in Sec. 2.3.2. Agreement with Eq. (2.30) entails < a < 0.6 
(compare also with Ref. [45], where less restrictive upper bound on w q (0) is imposed). The pa- 
rameter M can be determined for every a through Eq. (2.25) so that Eq. (2.29a) is satisfied. The 
determination of a and M is independent of ti, q\ and H\ provided that the tracking solution is 
reached in time. To reduce somehow the parameter space of our investigation we fix T\ - 10 9 GeV 
(or xi = -51.16). This choice is motivated by the majority of the inflationary models (see, e.g., 
Ref. [46]). We thereby focus on the two residual free parameters of our model and we design in 
Fig. 3-(a) [3-(b)] the allowed areas in the b - log Hi [b - log Q™ N ] plane. In the gray shaded areas 
Eqs. (2.26) - (2.30) are fulfilled. Obviously our model possesses an allowed parameter space with 
a band structure. The upper boundary curves of the allowed bands come from Eq. (2.28). Note, 
however, that saturation of Eq. (2.28) is not possible for 0.08 < b < 0.16. 

For any (b, log Hi), which is consistent with Eq. (2.28) and belongs in a white [gray] band, the 
resulting q after the oscillatory phase turns out to be negative [positive] and so, it cannot [can] serve 



attractor: 




(2.33) 
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as quintessence. E.g., let us fix b - 0.2. For 51.7 < \ogH\ < 53.3, q develops five extrema during 
its evolution - which is of the type shown in - resulting to qo > 0. Actually in Fig. 2-(b) we display 
the evolution of q as a function of x for the two bounds of this band. As log H\ decreases below 
53.3 (where the bound of Eq. (2.28) is saturated), the amplitude of the fifth peak, which appears 
in the ^-evolution (at about x - -24.5) eventually decreases and finally this peak disappears at 
\ogH\ - 51.7 where the first allowed band terminates. For 48.7 < \ogH\ < 51.7, q develops four 
extremes during its evolution resulting to q < 0. As log/Zj decreases below 51.7 the amplitude 
of the forth peak which appears in the ^-evolution (at about x - -30) decreases and finally this 
peak disappears at log Hi - 48.7 where the second allowed band commences. Note that, in the first 
allowed band, f2 BBN increases with Hi but this is not a generic rule. 

Variation of T\ (or equivalently ti) modifies the range of the obtainable f2 BBN and changes 
somehow the position of the several bands of our parameter space but do not alter the essential 
features of our results. E.g., for b = 0.2 and Ti = 10 10 GeV [Tj = 10 8 GeV] (or x Y = -53.5 
[ti = -48.9]) the margin of the first allowed band is 53.6 < log#i < 56.4 [49.6 < log#i < 51.1] 
and the second allowed band commences at log Hi - 50.7 [log Hi - 46.75]. 



3. The WIMP Relic Density 

We turn to the calculation of the relic density, Q. x h 2 , of a WIMP-CDM candidate, x- Employing 
the symbols defined in Eq. (2.9), Q. x h 2 can be found from the well-know formula: 



Q = PxO = ^o_ m x n x 

X PcO PcO S 



TYl f 

a x h 2 = 4.533 • 10~ 27 GeV 2 -^f 
x sR 3 



(3.1) 

Tf 



where s h 2 /p l c ^ = 4.533 • 10~ 27 GeV 2 and r f is a value oft [x] for the LRS [QKS] - see Eq. (2.10) 
[Eq. (2.21)] - large enough so as f x is stabilized to its present constant value, f x o. Recall that sR 3 
in the denominator of Eq. (3.1) remains constant only for x > trh in the LRS but for every x in 
the QKS - see Table 1. The evolution of f x obeys a Boltzmann equation. In Sec. 3.1 we present 
and solve this Boltzmann equation and in Sec. 3.2 we investigate how we can achieve within our 
scenaria an enhancement of Q. x h 2 w.r.t its value in the SC. 

3.1 The Boltzmann Equation 

Since ^'s are in kinetic equilibrium with the cosmic fluid, their number density, n x , evolves accord- 
ing to the Boltzmann equation: 

. ^„ ,(2 eq2\ ( ^<f N x n <P for the LRS, 

n x + 3Hn x + (av) (n x - n x ) = j Q ^ ^ QKS (3.2) 

where H is found from Eq. (2.1). In order to find a precise numerical solution to our problem, we 
have to solve Eq. (3.2) together with Eq. (2.20) [Eqs. (2.6a) and (2.6b)] for the QKS [LRS]. To this 
end, we rewrite Eq. (3.2) in terms of the quantities defined in Eq. (2.9) as follows 

HR 3 f' + 7^\(f 2 ^ 2 )- /^W? 6 for the LRS, 

HRf x + (crv)[f x -f x }-j for the QKS, (3 " 3) 

where H given by Eq. (2.12a) [Eq. (2.24)] and prime in this section denotes derivation w.r.t t = x 
[r = t] for the LRS [QKS]. Eq. (3.3) can be solved numerically in conjunction with Eq. (2.11) 
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[Eq. (2.23)] for the LRS [QKS]. In the LRS, we use the initial conditions in Eq. (2.13) and we 
integrate from x = to f f - 50 (with g's fixed to their values at Trh). In the QKS, we impose the 
initial condition h x {x x ) - tf x {x x ), where x x corresponds to the beginning (x = 1) of the Boltzmann 
suppression of n q in Eq. (2.7). The integration of Eq. (3.3) is realized from %\ down to Tbbn - 
-22.5 (an integration until to gives also the same result). 

Based on the semi-analytical expressions of Sec. 2, we can proceed to an approximate compu- 
tation, which facilitates the understanding of the problem and gives, in most cases, accurate results. 
In particular, we can set - see Eq. (2.24) [Eqs. (2.17) and (2.18)] for the QKS [LRS]: 

u F i, / 1 + r4 / r RH forx«XRH, , , I DC 

H- VPr^c where q c - < „ ,. for the LRS, or 

v [1 for x » trh 

1 [1+^Wpr forx«T K R, .... 

g c ~ \, H „ fortheQKS. (3.4) 

yL (1 - bq 2 /2) 1 1 for x » trr 

Introducing the notion of freeze-out temperature, 7p = T(j-p) - x?m x (see, e.g., Ref. [11, 16] 
and references therein) we are able to study Eq. (3.3) in the two extreme regimes: 

• For t <sc Tp, f x — f x q . So, it is more convenient to rewrite Eq. (3.3) in terms of the variable 
A(r) - f x (r) - J x \t) as follows: 

^(A- +r)+ ^A( i + 2^) = {W 6 ™ ; (3,, 

The freeze-out point tp can be defined by A(rp) = Spf^ q (Tp) where 5p is a constant of order 
unity, determined by comparing the exact numerical solution of Eq. (3.3) with the approximate 
under consideration one. Inserting this definition into Eq. (3.5), we obtain the equation: 

(* + l)^(r F) ^ (ln ^ (TF) + 6f{6f + 2) ^>^ (TF) = for the LRS, 

(i ^ / (16 + 3x)(18 + 25x) 

with (In/, ) ( T ) = 3 + x 2x2(g - 15x) (3-6) 

which can be solved w.r.t rp iteratively. The t - x [x - x] dependence can be derived by 
combining Eq. (2.7) and Eq. (2.2) [Eq. (2.12)] for the QKS [LRS]. 

• For t » t f , f x » f x q and so, we can set fjj: - J^ 2 - f£ in Eq. (3.3). Let us analyze this case 
for each scenario separately: 

♦ In the LRS and for the range of the parameters under consideration - see Sec. 3.2 - we 
single out two cases: 

* Dominant non-thermal production (non-TP). In this case, which is mainly applicable for 
very low Trh's, f x (crv) <s T^N^R 6 . Therefore, Eq. (3.3) can be integrated analytically 
inserting into it Eq. (2.17a) as follows: 

fxo - \ ^S^ x t^m-; (* 3W2 - e 3 ^ /2 ) . (3.7) 

Since f x o takes its main contribution close to trh » tp our result is more or less inde- 
pendent of fp. 
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* Equipartition between TP and non-TP Besides f x q none of the other terms of Eq. (3.3) 
can be neglected in this case. A rather precise result for f x o can be obtained by numeri- 
cally integrating Eq. (3.3) subsequently from (i) x = x p until x = xrh with initial condition 
f x {x ¥ ) = (<5 F + l)f x q {x ¥ ) and using H given by Eq. (3.4) [Eq. (2.17a)]; (ii) f - x RH un- 
til x = \f - 50, employing the estimates of Eq. (2.19) for and pr. A less accurate result 
for this case can be derived [9, 10] by equating the annihilation rate T x = n x {crv) to the 
expansion rate H at the completion of reheating. Combining Eq. (2.17a) and Eq. (2.18) 
we obtain //(x RH ) - 2p0(x RH ) - 2 V2f^/5 and so, we arrive at: 

Uy 9 V2 Trf, 

s n 2 (crv)Tl B 

In general, this result underestimates the numerical one by a factor of unity. However, 
the method applied reveals the presence of the phenomenon of reannihilation [47] in this 
case, i.e., the occurrence of a secondary (for x » xp)^ decoupling - see Sec. 3.2.1. 

♦ In the QKS, we can integrate numerically Eq. (3.3) from xp down to 0, as follows: 



fxo = (4"f + ^f)" 1 , where J F = £dxj^ and = (<5 F + 1) f^(x F ). (3.9) 



Although not crucial, a choice 5p = 1.2 + 0.2 assists us to approach better the precise 
numerical solution of Eq. (3.3). 

3.2 The Enhancement of Q. x h 2 

As we explain in Sec. 4 the interpretation of the e ± -CR anomalies favors 10~ 7 < (o~i>)/GeV~ 2 < 
10~ 6 which results to 0.0025 > & x h 2 \ sc > 0.00027 for 0.1 < rn^/TeV < 3, where & x h\ c denotes 
Q x h 2 within the SC. Clearly, the resulting Q x h 2 \ $c lie much lower than the range of Eq. (1.2). 
However, the proposed non-standard scenaria can increase Q x h 2 w.r.t Q^/i 2 | sc . The resulting en- 
hancement can be quantified, by defining the quantity: 

ao, = (a x h 2 - n x h 2 \ sc )/ a x h\ c . (3.10) 

We below analyze the behavior of AQ. X as a function of the free parameters of each non-standard 
scenario separately. 

3.2.1 The LRS. Let us initially clarify that in the LRS, both signs of AD. X are possible, as empha- 
sized in Ref. [11, 12]. However, we here confine ourselves to the combination of parameters which 
assures the favored from the e ± -CR data AQ. X > 0. The dependence of AQ^ on the free parame- 
ters of the LRS can be inferred from Fig. 4, where we depict AQ. X versus Trh for m x = 0.5 TeV, 
(o-v) = 10~ 6 GeV~ 2 [(av) = 10" 7 GeV~ 2 ] (gray [light gray] lines) and c x4> = 1 (solid lines), 
c x< p - 10~ 4 (dashed lines) and c x< p = 10" 6 (dotted lines). The ranges of parameters where each 
production mechanism is activated are also shown in the table included. Note that the exposed 
ranges depend very weakly on the employed m x % and (crvYs. We observe that AQ. X increases with 
Trh when we have non-TP as expected from Eq. (3.7), but it decreases as Trh increases when we 
have equipartition between non-TP and TP, as anticipated in Eq. (3.8). The former mechanism is 
dominant mainly for very low Trh's whereas the latter is present for higher Trh's. The accuracy 
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Figure 4: AQ. X versus Tm for the LRS with m x = 0.5 TeV, {crv) = 1(T 6 GeV~ 2 [((rv) = 1(T 7 GeV~ 2 ] (gray [light gray] 
lines) and c x< p = 1 (solid lines), c x $ = 10~ 4 (dashed lines) and c x< p = 10~ 6 (dotted lines). In the table we also show the 
type of thex production for each c x $ and the various ranges ofT m . 



of the corresponding empirical expressions in Eq. (3.7) [Eq. (3.8)] increases as Trh decreases [in- 
creases] and as (crv) decreases [increases]. It is remarkable that for c x $'s where both production 
mechanisms are possible (e.g., c x< p = 10~ 6 or 10~ 4 ) we can obtain the same AD. X for two values 
of 7rh- In general, AQ X increases with N x . Augmentation of m x increases £l x h 2 , too, but does 
not alter the dependence of Q. x h 2 on Jrh and the ranges where the ^-production mechanisms are 
activated. 

The operation of the two types of x production encountered in Fig. 4 is visualized in Fig. 5- 

(a) and (b). In these, we display the actual;^ yield, n x /s (solid lines) and its equilibrium value, 
n x js (dashed lines) the dimensionless actual interaction rate of x, = n x (crv) (solid lines), its 
equilibrium value, r^ q = ff x (crv) (dotted lines) and the dimensionless expansion rate H (dashed 
lines) - given by Eq. (2.12a) - versus f. In both figures we use m x = 0.5 TeV, - 100 TeV, 
(crv) = 3 • 10~ 7 GeV" 2 and (N X ,T RU ) = (1,0.5 GeV) [(N X ,T RU ) = (7.5 • 10 5 , 1 MeV)] (gray 
[light gray] lines). For the selected parameters, the evolution of the background energy densities 
(logp^ and logpR) and T is presented in Fig. l-(a) and (b). The completion of reheating occurs at 
T RH - 36.1 [T RH = 45.1] for T RR - 0.5 GeV [r RH = 1 MeV]. 

From Fig. 5-(a) we can deduce that n x /s takes its present value close to [clearly above] f RH for 
Trh = 1 MeV [Jrh = 0.5 GeV]. For this reason, the integration of Eq. (3.3) until trh for non-TP 
is sufficient for an accurate result - see Eq. (3.7) -, but insufficient when non-TP and TP interplay. 
The x reannihilation takes place along the almost vertical part of the gray line for f a little lower 
than f RH . It is notable that in Ref. [11, 12] which focus on lower (cru)'s than the ones considered 
here, the phenomenon of reannihilation is not stressed. 

This effect is further analyzed, following the approach of the first paper in Ref. [47], in Fig. 5- 

(b) . From this we infer that for Trh = 1 MeV, where non-TP outstrips, x decouples from plasma 
only once at Tf = 16.5 where T x - T x q - H and an intersection of the light gray lines is observed 
(note that the light gray and gray dashed lines coincide at that region). On the contrary, for Trh = 
0.5 GeV, where non-TP and TP coexist, we observe that^ decouples from plasma initially at about 
Tp = 26.3 where r^ q = H but also at Tp r = 37.9 > trh where T x = H. In other words, we observe 
two intersections between the two pairs of the three gray lines. This effect signalizes the existence 
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Figure 5: (a) n x /s (solid lines) and n^/s (dashed lines) versus x and (b) T x (solid lines), (dotted lines) and H (dashed 
lines) versus f for the LRS with m x = 0.5 TeV, m = 100 TeV, (crv) = 3 • 10~ 7 GeV~ 2 and (N x , T m ) = (1,0.5 GeV) 
[(N x ,T Rli ) = (7.5 ■ 10- 5 , 1 MeV)J (gray [light gray] lines) - in both cases, wegetSl x h 2 = 0.11. 



of a period of x reannihilation similar to that noticed in Ref. [47]. Contrary to that situation, in 
our case (i) T x remains larger than H after the first x decoupling and drops sharply below H after 
reannihilation, and (ii) H smoothly evolves from its behavior during LRS to that within RD era. 

3.2.2 The QKS. The presence of gc > 1 in Eq. (3.6) and, mainly, in Eq. (3.9) reduces /p w.r.t its 
value in the SC generating, thereby, AQ. X > within the QKS. The mechanism of the x decoupling 
in this case, for both b = and b + 0, is pretty known - see Ref. [8, 14, 16]. However, a peculiar 
effect emerges in the dependence of AQ^ on m x for b + which can be inferred from Fig. 6, where 
we display AQ^ versus m x for a = 0.5, H Y = 6.3 • 10 53 , (crv) = 1CT 6 GeV" 2 [(crv) = 1CT 7 GeV~ 2 ] 
(gray [light gray] lines) and b - (solid lines), b -0.15 (dashed lines) and b = 0.32 (dotted lines). 
The chosen Hi's result to Q^ BN 0.01, 0.068 or 0.19 for b = 0, 0.15 or 0.32 correspondingly. 

Obviously, for b = we get a pure KD era and our results reduce to those presented in 
Ref. [16], i.e., AQ. X increases when m x increases or (crv) decreases. On the contrary, for b + 0, 
b£l x depends crucially on the hierarchy between tp and T ex t found from Eqs. (3.6) and (2.32) 
respectively. Given that Jp takes its main contribution from gc for t ~ Tp, Jf is enhanced - see 
Eq. (3.9) - if tp is lower than x ex t and close to it, since gc is suppressed (gc - 1) for x - x ex t- As 
a consequence - see Eqs. (3.1) and (3.9) - AQ. X diminishes. This argument is highlighted in the 
table of Fig. 6. There, we list the range of xp for 0.1 < m^/TeV < 3 and (crv) = 10~ 7 GeV~ 2 or 
(crv) - 10~ 6 GeV and the logarithmic time x ex t at which the closest to xp's peak in the q evolution 
takes place for b = 0.15 or b - 0.32 and Hi = 6.3 • 10 53 . Clearly xp [x ext ] is independent of b and 
Hi [m x and (crv)]. As m x increases above 0.1 TeV, xp moves closer to x ext and A£l x decreases 
with its minimum AQ^| min occurring at xj? m - x ext . The small deviation of x™ n from x ext can be 
attributed to the presence of f X Y in Eq. (3.9). The appearance of the minima can be avoided if xp's 
happen to remain constantly lower than x ext 's - see, e.g., Fig. 9-(C2). Variation of Ti or Hi leads to a 
displacement of x ext 's - see Eq. (2.32) - relocating, thereby, the minima of ACl x in Fig. 6. However, 
our conclusions on the behavior of AQ. X remain intact. 
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Figure 6: AQ X versus m x for the QKS with a = 0.5, Hj = 6.3 • 10 53 , Ti = 10 9 GeV, (crv) = lO" 6 GeV~ 2 [{crv) = 
10~ 7 GeV -2 ] (gray [light gray] lines) and b = (solid lines), b = 0. 15 (dashed lines) and b = 0.32 (dotted lines). In the 
table we also show the range of x F for 0.1 < m^/TeV < 3, the freeze-out logarithmic time 1™" at which the minima of 
AQ. X occur and the closest to t™ ln 's x cxt 's for the selected b 's. 



4. PAMELA, ATIC and Fermi-LAT Anomalies 

The aforementioned AQ^ obtained within the LRS or QKS assists us to explain the experimental 
data on the e ± -CRs consistently with Eq. (1.2). Indeed, the observed anomalies on the CR e ± fluxes 
can be attributed to the annihilation of ^'s in the galaxy provided that m x and {crv) are chosen 
appropriately. In Sec. 4.1 we outline the basic formalism that we employ in order to estimate the 
observable quantities as a function of these parameters and in Sec. 4.2 we display our fittings. 

4.1 Cosmic Rays from Annihilation or WIMPs 

After being produced in the Milky Way halo, charged CRs propagate in the galaxy and its vicinity 
in a rather complicated way before reaching the earth. Their propagation is commonly evaluated by 
solving a diffusion equation [48, 49] with static cylindrical boundary conditions. The solution can 
be casted into the following semi-analytical form [48, 49, 50] which yields the e + flux per energy 
- in units GeV _1 cm^ 2 s _1 sr _1 - at earth from the^ annihilation: 

2Anb{E) \m x j J E dE' e+ 

where v e + is the velocity of e + which is practically equal to the one of the light, the pre-factor of 1/2 
arises from our assumption that^ is a Majorana particle, p Q = 0.3 GeV/cm 3 is the local CDM den- 
sity, b(E) - E 2 /(GeV ts) with tg - 10 16 s is the energy loss rate function and dN e +/dE e + denotes 
the energy distribution of <? + 's per^- annihilation. Motivated by the highly restrictive PAMELA 
data [51] on the anti-proton mode of the CDM annihilation, we consider a purely leptophilic x 
(however, see also Ref. [52]). In particular, we incorporate in our investigation the following an- 
nihilation modes: XX ~ > e e+ ^XX ~ > P or XX ~ > t~t + . In the first case, dN e +/dE e + is given 
by a simple analytic expression [53] whereas in both latter cases, we use analytic parametrizations 
(presented in Ref. [54, 55]) of dN e +/dE e + which reproduce quite accurately the numerical outputs 
of the package Pythia [56]. In all cases the effect of final state radiation [57] is taken into account. 
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Namely, we take 



dN e 
dE e 



for 
for 



XX 
XX 



(4.2) 



' 8{E e + - m x ) + (a sm /2n) ^36(E e + - m x )/2+ 
(l+y 2 )\n(4m//mfj/(l-y)m x ] 
(a em /n)A exp[-(A 1? / + A 2 y 2 )] + 5i + B 2 y 
(l/m^)(exp[A Q +A-y + A 2 y 2 +A-y 3 ] + 
, exp[A+ + A\y + A\y 2 + A+y 3 + A\y A + A+i/ 5 ]) for 

where a em is the fine-structure constant computed at a scale equal to 2m x 
e~ mass and < y - E e +/m x < 1. The infrared singularity encountered for XX — * e+e and (/ = 1 
is handled as described in Ref. [53]. For thexx — * mode, we take J - J(m x /0.5 TeV) dj with 
J = A, Ai, A 2 , Si, B 2 and 



AA ^ T T , 

m e = 0.511 MeV is the 



(A,A U A 2 ,B U B 2 ) = (-0.296635,2.65121, 14.8445,0.0042505,-0.00427157), 
(6 A , S Ai ,6 A2 ,6 Bi ,Sb 2 ) = (-1.01424, 0.017198, -0.0107585, -0.999819, -0.999819). 

We checked that the parametrization above gives results quite similar to those obtained using the 
simpler parametrization suggested in Ref. [55]. For thexx — * T+T mode, we take [55] 

(A , A[, A 2 , AJ) - (0.951, -18.083, 15.79, -15.575), 

(A^A^A^A^A^A^) = (2.783, -22.942, 82.595, -193.748,223.389, -97.716). 

Also, /(/Id) is the dimensionless halo function which fully encodes the galactic astrophysics 
with Ad(E, E') the diffusion length from energy E to energy E' which is given by 



l? D - 4K[)tE 



(E'/GeVy- 1 -(g/G eV) 



5-1 



(4.3) 



To compute I(ko) we employ the semi-empirical function proposed in Ref. [50, 58]. Namely, 

-I 



I(ko) - «o + a \ tanh 



ci 



, ( / ~^) 2 , , 
a 2 exp + a-$ 



with I - log ( 

\kpc 



(4.4) 



The involved in Eqs. (4.3) and (4.4) constants [50] depend on the CDM distribution and the prop- 
agation model that we consider. As we emphasize in Sec. 5.1.4, the constraint from the y CRs 
enforces us to adopt the isothermal halo profile [59] which weakens the relative restrictions. Note, 
however, that our results on 0*f are quite close to those that we would had obtained if we had 
used the NFW halo [60] profile - c.f. Ref. [31, 50]. We also use the MED propagation model for 
XX — > e e+ W &XX — > y u_ y u+ and the MIN (M2) model for XX ~ > t~t + . These choices provide the 
bets fits to the combined experimental data - c.f. Ref. [58, 61]. Note that only these two propaga- 
tion models are consistent [49] with the observed boron-to-carbon ratio in the CR flux. Therefore, 
we use [50] 



(ao,ai,a 2 ,a3,bi,b 2 ,ci,c 2 ,K Q ,6) - 
(0.495, 0.629, 0.137, 0.784, 0.766, 0.55, 0.193, 0.296, 0.01 12 kpc 2 /My, 0.7) (MED), 
(0.5, 0.903, -0.449, 0.557, 0.096, 192.8, 0.210, 33.91, 0.00595 kpc 2 /My, 0.55) (MIN). 



(4.5) 



We explicitly verified that the numerically fitted function in Eq. (4.4) reproduces quite accurately 
and fast enough the results obtained by performing numerically the relevant integrations presented 
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in the earlier formulae of Ref. [49]. Moreover, the formalism of Ref. [50] overcomes successfully 
the mismatching problem (in the numerical integration) which is mentioned in Ref. [62]. 

In order to calculate the total fluxes, we also have to estimate the background e ± fluxes. In our 
study, we take into account the fluxes of (i) secondary e + (Q>^f) produced by collisions between 
primary protons and interstellar medium in our galaxy (ii) primary e~ (o^ nm ) presumably produced 
in supernova remnants and (iii) secondary e~ (^g 20 ) produced by spallation of CRs in the interstellar 
medium. These fluxes are commonly parameterized as [48] 

sec = 4.5(£/GeV)°- 7 

1 + 650 (£/GeV) 2 - 3 + 1500 (£/GeV) 4 - 2 ' 

<D*» = 0.16(g/GeVT" 

1 + 11 (E/GeV) - 9 + 3.2(£/GeV) 2 - 15 ' 

1 + 1 10 (E/GeV) 1 - 5 + 600 (£/GeV) 2 - 9 + 580 (£/GeV) 4 - 2 ' 
in units GeV _1 cnr 2 s _1 sr _1 . With these backgrounds, the total e ± fluxes read 

dV = ®£ + 0>f + c and O e - = + c e - <$> v ™ + 0* ec (4.7) 

where & x = and c e - (0.6 - 0.8) is a normalization factor. We take c e - = 0.6 [c e - - 0.7] so 
that our flux calculation is consistent with the ATIC [Fermi-LAT] data in the low energy range of 
(20 - 70) GeV [58]. 

4.2 Fitting the PAMELA and ATIC or Fermi-LAT Data 

Using the fluxes defined above, we can evaluate the observable quantities and compare them with 
the experimental outputs. In order to qualify our fittings to the experimental data, we perform a % 2 
analysis. In particular, we define the % 2 variables as [54, 58, 62] 

( F obs _ F £\ 2 1 0g+ / ( o e+ + O e ) and N A = 7 for A = 1, 

2 \ v_ai AO F \ (21 for A -2, (4.8) 

(AFff for A -_ x 



i=l 



where A = 1,2,3 stands for the PAMELA [24], ATIC [25] or Fermi-LAT [30] data respectively 
which are considered as independent sets. The index i runs over the data points of each experiment 
A, the superscript "obs" ["th"] refers to measured [theoretically predicted] quantities whereas AF obs 
means error in the experimentally observable F. Na is the number of data points considered from 
the experiment A. Note that, from the PAMELA data-set, we use [58, 54] only the 7 data points 
above 9.1 GeV where the effect of solar modulation is expected to be small. In our analysis we 
take into account only the vertical errors. We also conservatively combine, independently for each 
data-point, in quadrature statistical and systematic errors released from Fermi LAT [30]. 

In Fig. 7 we show the predicted observable quantities compared to the experimental data as a 
function of the e + energy E e +, assuming x annihilating to e + e~ (dot-dashed lines), fi + /u~ (dashed 
lines) or t + t~ (dotted lines). We use the best-fit (m x , (<ru))'s obtained from minimization of % 2 + yjj, 
[X? + X3] m Fig- 7-(a0 and {a.2) [Fig. 7-(bi) and (D2)]. Since ATIC and Fermi-LAT data are not 
consistent with each other, we do not combine them but present results using only either of the two. 
The relevant (m x , (cry»'s can be read in the Table 2 together with the corresponding x 2 /d.o.f, where 
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Figure 7: The e + -3ux fraction (ai and bj and the total e~ and e + flux times E? + (a. 2 and t> 2 ) as a function of E e +, with 
E e + being the e + energy. We use the best-fit points (m x , (crv)), indicated in Table 2, for the PAMELA and ATIC data (ai 
and ai) or the PAMELA and Fermi-LAT data (bi and b 2 ), assuming x annihilating into e + e~ (dot-dashed lines), \i*pT 
(dashed lines) orr + T~ (dotted lines). The background fluxes are denoted by solid lines and are computed for c e - = 0.6 
(ai and ^2) or c e - = 0.7 (bi and b2). The data from PAMELA (ai , bi), ATIC (82) and Fermi-LAT (b2> experiments are 
also shown (an additional uncertainty from the Fermi-LAT energy scale, which can shift all the points by 5% (up) to 
10% (down) is not shown). Recall that we adopt the isothermal halo profile and the MED [MIN] propagation model for 
XX -» e + e- and xx -» H*H~ kx -» 

d.o.f denotes the number of degrees of freedom involved in our fits which is equal to N1+N2-2 - 26 
[Ny + N 3 - 2 - 31] for PAMELA and ATIC [PAMELA and Fermi-LAT] data (2 is the number of 
the fitting variables, m x and {crv)). As can be deduced from Table 2, an exceptionally good fit to 
PAMELA and Fermi-LAT data arises for^ annihilating to yu + ju _ , whereas in most other cases the 
fits are rather poor since we get just x 2 /d.o.f ^2.5-3 for 28 or 33 data points. Better fits can 
be probably attained under the assumption that ^'s both annihilate and decay, as pointed out in 
Ref. [58]. Note finally that, the^ annihilation into e + e~ is strongly disfavored [63] by Fermi-LAT 
data since it predicts a spectrum with a too sharp end-point. 

In Table 2 we also list the maximal {crv), {crv) max , allowed by Eq. (5.4) (see Sec. 5.1.3) and 
the resulting 0. x h 2 \ 's. It is remarkable that O^/i 2 ] turns out to be well below the range of 
Eq. (1.2) implied by the CDM considerations. As a consequence, the SC can not be consistent with 
the interpretation of the e ± -CR anomalies via x annihilation, unless we invoke an enhancement 
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Mode 
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Q.O.I 


lev 


1U ueV 


1U uev 


1 a -4 
1U 


7-(a0 and (a 2 ) 


XX -» 


67/26 


0.74 


7 


3.3 


3.85 


7-(ai) and (a 2 ) 


XX -» 


76/26 


2 


28.6 


26 


0.97 


7-(a0 and (a 2 ) 


XX -* T + T- 


79/26 


3.5 


143 


20 


0.2 


7-(bi) and (b 2 ) 


XX -* e + e~ 


75/31 


0.5595 


4.6 


2.5 


5.8 


7-(bi) and (b 2 ) 


XX -» 


24/31 


1.28 


19.5 


16.5 


1.4 


7-(a0 and (a 2 ) 


A'A' -» 


33/31 


2.1 


94.6 


28.2 


0.3 



Table 2: The best-fit points (m x , (crv)) for the PAMELA and ATIC data (used in Fig. 7-(a,) and (a 2 )) or the PAMELA 
and Fermi-LAT data (used in Fig. 7-(t>i ) and (b2> ) for each annihilation channel ofx- Reported are also the corresponding 
X 2 | min /d.o.f ., the maximal (crv), (<ru} raax , allowed by Eq. (5.4) and the resulting £l x h 2 in the SC, Q. x h 2 \ SQ . 



mechanism of (crv) at present [26, 27]. In other words, some AQ. X is necessitated in order to 
reconcile the best-fit (m x , (<xu))'s with Eq. (1.2). Moreover, we observe that the bound of Eq. (5.4) 
(which turns out to be the most restrictive of all the others presented in Sec. 5.1) is violated in 
all cases. This violation is softer [stronger] in the case where ^'s annihilate to /i + ju [t + t~]. We 
observe, also, that for XX — * T+T ~> "V s an d (cn>)'s are pushed to larger values than those needed 
f° r XX ~ > A^j" an d the m x \ and (cn;)'s used in the latter case are higher than those used for 
XX — > e + e~. As shown in the next section, best-fit (m x , (<xu))'s consistent with all the available 
constraints can be achieved fork's annihilating to e + e~ or /i + ju~. 

5. Restrictions in the m x - (crv) Plane 

To systematize our approach, we need to delineate in the m x - (crv) plane the regions which are 
favored at 95% c.l. by the various experimental data on the e ± -CRs. In particular, we consider 
[58, 54] regions favored by PAMELA data only, PAMELA and ATIC data or PAMELA and Fermi 
LAT data. Since we have two independent parameters, m x and (crv) these regions can be determined 
imposing the condition [54] 

(- Xi for PAMELA, 

X 2 * X 2 L + 6 with X 2 = Xi + X2 for PAMELA and ATIC, (5.1) 

V X 2 + X3 for PAMELA and Fermi LAT, 

where X 2 | min can De extracted numerically by minimization of % 2 w.r.t m x and (crv). On the other 
hand, the interpretation of the data on e ± -CRs in terms of x annihilation can be viable if it can 
become consistent with a number of phenomelogical constraints. In Sec. 5.1 we summarize these 
constraints and in Sec. 5.2 we examine if they can be reconciled with the regions favored by the 
data on e ± -CRs. 

5.1 Imposed Constraints 

Though n x /s in Eq. (3.1) stays essentially unchanged for r > Tf, residual annihilations of _^'s 
occur up to the present with several cosmological consequences besides the possible interpretation 
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of the data on e ± -CRs. Recently, important upper bounds on (crv) have been reported and are 
summarized below for the three exemplary x annihilation modes considered in our investigation. 
The well-known unitarity constraint is also taken into account. 

5.1 .1 Unitarity Constraint. Using partial-wave unitarity [23, 65] an upper limit, particularly rel- 
evant for m x > 2 TeV, on (crv) can be derived as a function of m x , i.e., 

(in \ ~2 
IGev) • (5 " 2) 

5.1 .2 BBN Constraint. During BBN, the^- annihilations inject an amount of energetic particles 
which is proportional to (crv) and may strongly alter [61, 64] the abundances of the light elements. 
Ruining the successful predictions of the BBN can be avoided if we impose an upper bound on 
(crv) which, however, depends on the identity of the products of the annihilation of _^'s. Taking into 
account the most up-to-date analysis of Ref. [61] we demand: 

2m m E f 2 for XX ~» « + «", 

(crv) < 3 • 1(T 5 GeV" 2 -r^—^ where = 0.7 for X X -» ft, (53) 
£ vis 1 TeV m x 

A \ 0.62 for XX T T ' 

with £" v i s being the total visible energy of the produced particles in the X annihilation. 

5.1 .3 CMB Constraint. The^- annihilations may have [32, 33] an impact on the ionization state 
of the baryonic gas at recombination and therefore on the CMB angular spectra. Consistency with 
the WMAP5 data [1] dictates [33] at 95% c.l. (see, also, Ref. [34]): 

-7 o f0.7 for vv^>e + e~, 

3.1 • 10" 7 Ge^ 2 m y AA 
(crv) < -f- where / - 0.24 for X X -» y"V, (5-4) 

1 ' ' 1 0.23 for xx^r + r-, 

is the deposited power fraction which expresses the efficiency of the coupling between the annihila- 
tion products and the photon-baryon fluid at z ~ 1000. It is expected that forthcoming experiments 
will impose [33] even more stringent bounds on (crv). Note, in passing, that the presence of q in the 
QKS does not affect recombination (which occurs at x rec - -7) since £l q (x mc ) is safely suppressed 
provided that Eq. (2.28) is fulfilled. 

5.1 .4 Constraint from the y-CosMic Rays. The^f annihilation in the galactic center yields sizeable 
amount of y-CRs, through the cascade decay of the annihilation products and/or bremsstrahlung 
processes. Comparing the relevant y-CR flux with the H.E.S.S observations [66] we can further 
restrict [31] (crv) as a function of m x for the two chosen annihilation channels. However, this 
restriction significantly depends on the CDM halo profile. Adopting the cored isothermal CDM 
profile [59], which assures the less restrictive version of this constraint, we graphically extract the 
upper bound on (crv) from the plots of Ref. [31]. To have a feeling of the strength of this constraint 
we can give some rough estimations: 

(- (3 - 10) • 10" 6 GeV" 2 for m x - (0.2 - 1.5) TeV and X X -» e + e", 
(crv) < | (4 - 25) • 10~ 6 GeV" 2 for m x - (0.2 - 3) TeV and X X -> A<V, (5.5) 

{ (1.2 - 4) • 10" 6 GeV" 2 for m x = (0.5 - 3) TeV and X X -> t + t~- 
Alternatively this constraint can be evaded for every CDM profile, by allowing the ^-annihilation 
products to be long lived, as pointed out in Ref. [54]. 
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Complementary constraints on the x annihilation can be imposed comparing the findings of 
EGRET satellite [67] with the diffuse (secondary) y-CR fluxes, which would be produced [68] by 
inverse Compton scatterings on interstellar photons of the energetic e ± generated by the x annihi- 
lation in the galactic halo. However, these bounds are expected [68] to be weaker than the ones 
imposed by the high energy y-CRs mentioned in Sec. 5.1.4 and are not included in our analysis. 
Similar arguments are [69] also valid for the neutrinos generated from the x annihilation in the 
galactic center, though the dependence on the CDM profile is weaker. 

5.2 Results 

Constructing the preferred areas by the various experimental data and taking into account the con- 
straints quoted in Sec. 5.1 we can check the viability of the interpretation of the anomalies on 
e ± -CR fluxes in terms of the^f annihilation. In Fig. 8 we consider the mode^ — > t + t _ whereas in 
Fig. 9-(a0, (bi) and (ci) [Fig. 9-(a2), (b2) and (C2)] we assume mate's annihilate into e + e~ [p^^T]. 

In Fig. 8 and 9 we delineate the regions preferred at 95% c.l. by PAMELA data (black and 
red sparse hatched areas), PAMELA and ATIC data (dense black hatched areas) and PAMELA 
and Fermi-LAT data (dense red hatched areas), by imposing the condition of Eq. (5.1). In the 
black [red] hatched areas the backgrounds fluxes are normalized setting c e - - 0.6 [c e - = 0.7] in 
Eq. (4.7). Evidently, the PAMELA data do not prefer any m x since it does not show any peak 
structure whereas Fermi-LAT data disfavors the mode XX ~ * e+e since the spectrum from such 
a channel is too peaked to reproduce data. We also remark that the regions derived by the joint 
analysis of two data-sets are rather limited - c.f. Ref. [31, 34, 54, 63]. We consider the latter results 
as more reliable, since even when we fit only the PAMELA data, the data points with low E e + from 
ATIC/Fermi LAT are involved, in order to normalize [54] the background fluxes. 

In Fig. 8 and 9 drawn is also the upper bound from Eqs. (5.2), (5.3) and (5.4), denoted by a 
black dotted, solid and dashed line respectively and this from the constraint of Sec. 5.1.4, depicted 
by a dot-dashed line. Obviously acceptable are the regions mainly below the dashed curves, since 
the bound of Eq. (5.4) is the most restrictive from the others. The bound of Eq. (5.2) cuts out some 
slices of the parameter space fov xx — * an ^XX — * T+T an d large m x s. We easily conclude 
that the explanation of the experimental anomalies via the annihilation mode: 

• XX e+e ~ is just marginally consistent with Eq. (5.4) - see Fig. 9-(a0, (bi) and (ci). Indeed, 
we observe that just a minor portion of the area favored by PAMELA at 95% c.l. is allowed by 
Eq. (5.4) whereas the regions preferred at 95% c.l. from both combinations of PAMELA and 
ATIC or PAMELA and Fermi-LAT data are entirely excluded from the bounds of Eq. (5.4). 

• XX can be reconciled - c.f. Ref. [23, 31, 54, 61] - with the various constraints 
- see Fig. 9-(a2), (b2) and (C2). Namely, we notice that sizable slices of the regions favored 
by PAMELA lie below the bound of Eq. (5.4). Moreover, very close to or even lower than 
this limit we find portions of the favored regions at 95% c.l. by the PAMELA and ATIC or 
Fermi-LAT data. 

• XX ~ > T+T ~ is inconsistent - c.f. Ref. [31] - with both Eq. (5.4) and Eq. (5.5) at 95% c.l. 
since all the regions favored by the experimental data lie entirely above the bounds above - 
see Fig. 8. Violation of the bound of Eq. (5.2) in a sizable fraction of these regions is observed 
too. Because of this fact, we below concentrate on the other two annihilation modes of x's- 
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Figure 8: Restrictions in the m x - (crv) plane 
forx's annihilating into t + t~ . The sparse black 



> 

1 r 



XX- T T 



[red] hatched areas are preferred at 95% c.l. by 
the PAMELA data for c e - = 0.6 [c e - = 0.7J 
and the dense black [red] hatched areas are 
preferred at 95% c.l. by the PAMELA and 
ATIC [PAMELA and Fermi-LAT] data. Re- 
gions above the black solid, dashed, dotted and 
dot-dashed lines are ruled out by the upper 
bounds on (crv) from Eq. (5.3), (5.4), (5.2) and 
Sec. 5.1.4 correspondingly. 
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Fixing the parameters related to the LRS [QKS] to some representative values - consistent 
with Sec. 2.2.2 [Sec. 2.3.2] -, we can display in the m x - (crv) plane, as in Fig. 9, regions (light 
gray shaded) confronted with Eq. (1.2). The gray dashed [dotted] lines correspond to Eq. (1.2b) 
[Eq. (1.2a)], whereas the gray solid lines are obtained by fixing Q. x h 2 to its central value in Eq. (1.2). 
Fig. 9-(a0 and (a 2 ) are devoted to the LRS whereas Fig. 9-(t>i), (b 2 ), (cO and (C2) analyze the QKS. 

For the LRS, we confine ourselves to some combinations of parameters which assure a suffi- 
cient coexistence of non-TP and TP, since non-TP alone is obviously - see Eq. (3.7) - (crv) inde- 
pendent and therefore can not be properly depicted in the m x - (crv) plane. We take (c x $, 7rh) = 
(2 • 10" 6 ,0.1 GeV) and (c^,r RH ) - (1, 1 GeV) [(c^,r RH ) - (10" 6 ,0.1 GeV) and (c^,r RH ) - 
(1,0.5 GeV)] in Fig. 9-(ai) [Fig. 9-(a 2 )]. For the QKS, we set throughout a = 0.5, Tj = 10 9 GeV. 
We also take (i) b = and Hi = 6.3 • 10 53 , 2 • 10 53 or 6.2 • 10 52 resulting to H BBN = 0.01,0.001 
or 0.0001 respectively in Fig. 9-(bi) and (b 2 ); (ii) Hi = 6.3 • 10 53 and b = 0.15 [Hi = 6.2 • 10 52 and 
b = 0.32] resulting to Q BBN = 0.068 or [Q BBN = 0.065] in Fig. 9-(Ci); (iii) Hi = 6.3 • 10 53 and 
b = 0.32 [#1 = 2- 10 53 and b = 0.2] yielding fi BBN = 0.19 or [Q BBN - 0.21] in Fig. 9-(c 2 ). Note 
that in Fig. 9-(bi) and (b 2 ) we present for the sake of comparison results even for b = 0, although 
the tracking behavior of the QKS is not attained in this case - see Sec. 2.3.4. 

In all cases, we observe that Cl x h 2 decreases as (crv) increases. This is due to the fact that 
n x h 2 oc 1 1 (crv) as can be deduced from Eq. (3.1) and Eq. (3.8) [Eq. (3.9)] for the LRS [QKS]. 
For the LRS , as it is clear from these plots, there is a minor slice of the allowed region with 
m x < 0.35 TeV [m x < 0.8 TeV ] for r RH = 0.1 GeV and c x<p = 2 ■ 10~ 6 [c x<p = 10~ 6 ] where non-TP 
is strengthened and our results are almost (crv) independent. For the QKS, we also observe that 
for tp far away from x ext the allowed by Eq. (1.2) for b + region reaches the one for b = - 
with fixed H\. However, when tp reaches x ext , Q. x h 2 decreases (as we explain in Sec. 3.2.2) and so, 
the required, for obtaining Cl x h 2 in the range of Eq. (1.2), (crv) decreases too. As a consequence, 
although the allowed by Eq. (1.2) area in Fig. 9-(ci) [Fig. 9-(c 2 )] for Hi = 6.3 • 10 53 approaches the 
corresponding area in Fig. 9-(bi) [Fig. 9-(b 2 )] with the same Hi and violates the bounds of Eq. (5.4) 
for low m x 's, it becomes compatible with the latter constraint for larger m x s. On the other hand, we 
observe that there is no such a transition region in the light gray area of Fig. 9-(c 2 ). This is, because 
for 0.1 < ra^/TeV < 3 we get 31.8 < -Xp < 35.1 whereas the closest to xp's, x ext is x ext = -35.4 
which remains constantly lower than xp. Therefore, no reduction of A£l x occurs for the m x s used 
in Fig. 9-(c 2 ). 
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Figure 9: Restrictions in the m x - (crv) plane for the LRS (ai, a2> [QKS (b 



taking a = 0.5, T\ = 10 GeVJ with various c x $ 's and Trh 's [b's and Hi 's], indicated in 
the graphs, andx's annihilating into e + e~~ (ai, bi, Ci) or (ai, b2, C2). The light gray 
shaded areas are allowed by Eq. (1.2), the sparse black [red] hatched areas are preferred at 
95% c.J. by the PAMELA data forc e - = 0.6 [c e - = 0.7J and the dense black [red] hatched 
areas are preferred at 95% c.J. by the PAMELA and ATIC [PAMELA and Fermi-LAT] 
data. Regions above the black solid, dashed, dotted and dot-dashed lines are ruled out by 
the upper bounds on (crv) from Eq. (5.3), (5.4), (5.2) and Sec. 5.1.4 correspondingly. The 
conventions adopted for the residual lines are also shown. 
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As can be concluded from most of the plots of Fig. 9, a simultaneous interpretation of the 
e ± -CR anomalies consistently with the requirements of Sec. 5.1 can be achieved in the regions 
where the gray shaded areas overlap the lined ones below the dashed lines. To clarify further this 
intriguing conclusion of this paper, it would be interesting to find the best-fit (m x , (cry)) (for the 
various combined data-sets) which fulfill all the restrictions imposed in Sec. 5.1. Our results are 
arranged in Table 3. The listed (m x , (o~u))'s saturate the bound of Eq. (5.4) which turns out to 
be essentially the most stringent of the others - see Fig. 9. We observe that m x ~ 0.1 TeV and 
{crv) ~ 1(T 7 GeV 2 [m x ~ 1 TeV and (crv) ~ 1(T 6 GeV~ 2 ] fork's annihilating to e+e' \ji + fi~l 
From the exposed in Table 3 x 2 - X 2 | m i n ' s » we deduce that all the requirements are met in a portion 
of the area favored at 99% c.l. [68% c.L] by the PAMELA and Fermi-LAT [ATIC] data for 
XX whereas the modexx — * e+e ~ can be excluded at 99% c.l. As regards the quality of 

the fits, from Tables 2 and 3, we can infer that the p^yT channel gives better fit to the PAMELA 
and Fermi-LAT data (x 2 /d.o.f = 33/31) than to the PAMELA and ATIC data 0< 2 /d.o.f - 77/26). 

From Table 3 we can also appreciate the importance of the non-SC in boosting £l x h 2 to an 
acceptable level. Indeed, in this Table we display 0^/i 2 | sc for every allowed best-fit (m x , (crv)). We 
observe that it lies much lower than the range of Eq. (1.2) in all cases, i.e., it is insufficient to account 
for the present CDM abundance in the universe. However, an appropriate adjustment (shown also 
in Table 3) of the parameters c X( p and Trh [b and Hi (with fixed a = 0.5 and T\ = 10 9 GeV)] for 
the LRS [QKS] - consistently with the restrictions of Sec. 2.2.2 [Sec. 2.3.2] - elevates adequately 
Q. x h 2 which can acquire the central experimental value in Eq. (1.2) for every best-fit point. In 
Table 3 we also expose the type of x production for the LRS and the transition temperature to the 
RD era for the QKS. We remark that since the (cr)'s required fork's annihilating to e + e~ are lower 
than those required for the y^yT channel, non-TP dominates even for Trh = 0.1 GeV. Note that 
Tkr < 0.04 GeV and the tracking behavior fails for b - in the case of the QKS. 

6. Conclusions 

We presented two non-standard cosmological scenaria which can increase the relic abundance of 
a WIMP x, Q x h 2 , w.r.t its value in the SC due to the generation of a background energy density 
steeper than the one of RD era. This increase is quantified by A£l x defined in Eq. (3.10). According 
to the first scenario, termed LRS, a scalar field (f> decays, reheating the universe to a reheating tem- 
perature lower than the freeze-out temperature of the WIMPs. According to the second scenario, 
termed QKS, a scalar field, q, rolls down its inverse power-low potential with a Hubble-induced 
mass term. In both cases our approach was both (i) purely numerical, integrating the relevant system 
of the differential equations (ii) semi-analytical, producing approximate relations for the evolution 
of the various energy densities of the cosmological background and the ^-number density. We 
consider that the exposed semi-analytical findings - although do not provide quite accurate results 
in all cases - facilitate the understanding of the cosmological dynamics. 

As regards the LRS, we recalled the dynamics of reheating and showed that AQ. X is affected 
by the two basic types of x production which can be discriminated, depending whether non-TP 
dominates or equally contributes with TP. The first type is activated for very low Trh, low N x and 
is more or less independent of (crv), whereas the latter case requires larger 7rh's and N x s and 
dependents on (crv). In this last case, we remarked that a period of x reannihilation can emerge. 
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Fits to PAMELA And ATIC Data 
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I 1 -I 2 


Imin 


13.4 


1 
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Combinations of Parameters Yielding Q. x h 2 


= 0.11 In The LRS 
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0.001 
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Fits to PAMELA And Fermi-LAT Data 
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Combinations of Parameters Yielding Q. x h 2 


= 0.11 In The LRS 




7r H (GeV) 


0.001 


0.1 
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0.1 
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Combinations of Parameters Yielding Cl x h 2 =0.11 








In The QKS for a = 0.5 and T\ = 


10 9 GeV 






b 







0.2 


0.32 





0.08 


0.18 


H^IO 53 


0.79 


1.3 


0.99 


3.1 


3.4 


4.7 


Tkr (GeV) 


0.07 


0.03 


0.04 


0.019 


0.005 


0.009 



Table 3: Best-fit (m x , {crv))'s for the combination of the PAMELA and ATIC or Fermi-LAT data and the various 
annihilation channels, consistently with all the imposed constraints. Shown are also the resulting / 2 - X 2 | min an d Q x h 2 
in the SC, Cl x h 2 \ , several combinations of(b,H{)'s [(T m ,c x4 ,)'s] leading to Cl x h 2 ^0.11 and the corresponding Tkr's 
[types of x production] within the QKS [LRS]. 
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As regards the QKS, we verified that the included Hubble-induced mass term ensures the 
presence of a KD period, which is characterized by an oscillating evolution of q, and allows the 
quintessential energy density to join in time a tracker behavior, alleviating, thereby, the coincidence 
problem. Observational data originating from BBN, the present acceleration of the universe, the 
inflationary scale and the DE density parameter can be also met in a sizable fraction of the parame- 
ter space of the model. AQ. X crucially depends on the hierarchy between the freeze-out temperature 
and the temperature where the evolution of q develops extrema. 

Assuming that the WIMP annihilates primarily to e + e~, ju + ju~ or t + t~ we calculated the in- 
duced flux of f^-CRs and fit the current data of PAMELA, ATIC and Fermi LAT without invoking 
any ad-hoc boost factor. For simplicity, we did not include in our fits older experimental results, 
such as from PPB-BETS [70], or more uncertain ones, such as from H.E.S.S [71] (however, the 
latter data may be used for imposing an upper limit in the m x - (crv) plane [72]). Taking into ac- 
count the strong bounds originating mostly from CMB, we concluded that the channel^ — > t + t~ 
can be excluded at 95% c.l. and: (i) large parts of the regions favored by PAMELA at 95% c.l. for 
the residual annihilation modes are ruled out; (ii) regions favored by PAMELA and ATIC or Fermi 
LAT at 99% c.l. for XX ~ * e+e 316 excluded; (iii) only a part of the region favored by PAMELA 
and ATIC data at 95% c.l. forxx ~ * ^ + ^~ can be acceptable. For the latter annihilation channel we 
achieved our best fits to PAMELA and Fermi-LAT data with m x ~ 1 TeV and (crv) ~ 1(T 6 GeV -2 
which belong within the region individuated by these data-sets at 99% c.l. In all cases, the consid- 
ered m x s and (cru)'s can yield the right amount of Q. x h 2 (entailed by the CDM considerations) by 
adjusting the parameters of the QKS or LRS. In both scenaria the required transition temperature 
to the conventional RD era turns out to be lower than about 0.7 GeV. In the case of the LRS an 
appreciable contribution of non-TP is also necessitated. 

As for the prediction of any CDM signal, there are three sources of uncertainty in our in- 
vestigation: the CDM distribution, the propagation of its annihilation products and the role of 
astrophysical backgrounds. In our analysis we used (i) the isothermal halo profile, to avoid troubles 
[31] with observations on y-CRs; (ii) the MED propagation model, which provides the best fits to 
the combinations of the various data-sets [58] and (iii) commonly assumed background <? + and e~ 
fluxes [48] normalized with the ATIC or Fermi-LAT data. The uncertainties above in conjunction 
with the very stringent constraints from CMB [33, 34] may jeopardize the interpretation of the 
PAMELA and ATIC or Fermi-LAT anomalies through the x annihilation. Therefore, the proposed 
scenaria can be probed in the near future, if a better understanding of the astrophysical uncertainties 
becomes available and/or more accurate experimental data are released. 

Our proposal could be supplemented by the construction of a particle model (see, e.g., Ref. [58, 
73]) with the appropriate couplings so that^f annihilates into ju + /i~ with the desired (trie's derived 
self-consistently with the (s)particle spectrum. In a such case, several phenomenological implica- 
tions could be examined as in Ref. [15]. Let us finally mention that another class of non-standard 
cosmological scenaria [8, 23, 47, 74, 75] can be generated considering modifications to the Fried- 
mann equation due to corrections to the Einstein gravity. Constraining these possibilities (as in 
Ref. [76]) through the experimental results on e ± -CRs would be another interesting issue. 
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